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SUMMARY 


An effort to develop a three-dimensional ice accretion modeling method has been 
initiated. This first step towards creation of a complete aircraft icing simulation code 
builds on previously developed methods for calculating three-dimensional flowfields and 
particle trajectories combined with a two-dimensional ice accretion calculation along 
coordinate locations corresponding to streamlines. This work is intended as a 
demonstration of the types of calculations necessary to predict a three-dimensional ice 
accretion. Results of calculations using the 3D method for a MS-317 swept wing 
geometry are projected onto a 2D plane parallel to the free stream direction and 
compared to experimental results for the same geometry. It is anticipated that many 
modifications will be made to this approach, however this effort will lay the groundwork 
for future modeling efforts. Results indicate that rime ice shapes predicted by the code 
agree very well with measured ice shape profiles. Predictions of glaze ice shapes 
indicate a difficulty in accurately calculating the ice shape in the runback region. 


NOMENCLATURE 


A m 

= area that particles impact on the 

h, HTC = 

convective heat transfer 


surface 


coefficient 

A 0 

= area through which a group of 

i. j = 

grid indices 


particles travel at release plane 

K 

latent heat of vaporization 

A s 

= control volume surface area 

m c 

mass flux determined from 

C P 

= pressure coefficient 


collection efficiency values 

C P 

r w 

= specific heat of water 

m e 

mass flux evaporating from the 

’Lee 

= ice thickness 


surface 

/, FF 

= freezing fraction 

mj 

mass flux freezing 
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m r . 

1 in 

= 

runback mass flux from the 

T s , TS 

~ 

temperature at the surface 



previous control volume 

T oo 

= 

temperature in the free stream 

m r 

*out 

= 

mass flux into the downstream 

t 

= 

time 



control volume 

u 

= 

weighting factor in j-direction 

<ic 

= 

heat flux due to the impacting 

VE 

= 

local inviscid velocity 



water droplets 

Voo 

= 

free-stream velocity 

^cond 

= 

heat flux due to conduction 

V 

= 

weighting factor in i-direction 



through the ice/airfoil surface 

w 

= 

unit width in spanwise direction 

Qconv 

= 

heat flux due to convection 


= 

displacement vector from x(i j) to 

Qevap 

= 

heat flux due to evaporation 



x(i+lj) 


= 

heat flux required for phase 


= 

displacement vector from x(i j) to 



change of water to ice 


x(ij+l) 

**in 

= 

heat flux associated with the 


= 

displacement vector from x(i j) to 



incoming runback water 



x . 

^ out 

= 

heat flux associated with the 

x . 

= 

location along surface streamline 



outgoing runback water 

x c 

= 

location of centroid of A m 

Re 

= 

Reynolds number 

/?,BETA 

= 

collection efficiency 

S 

= 

surface 

7 

= 

sweep angle 

s 


surface distance 

P 

= 

density 

S C 

= 

chordwise component of s 





I. Introduction 

Ice accretion modeling has received considerable attention during the past few 
years. The National Aircraft Icing Technology Plan 1 , drafted in 1986 under the 
direction of several federal agencies, calls for the creation of a complete aircraft icing 
simulation during the coming decade. In support of this plan several activities related 
to the creation of a 2D ice accretion modeling and performance evaluation capability 
have been undertaken at the NASA Lewis Research Center. A 2D ice accretion code, 
LEWICE, was developed in 1983 by the University of Dayton Research Institute 2 and 
later modified by Ruff 3 . This code is divided into three major parts; a potential flow 
calculation, a particle trajectory calculation, and an ice growth calculation. The ability 
of LEWICE to predict airfoil ice accretions is quite good for rime ice growths. Glaze ice 
predictions, although acceptable, have exhibited some deficiencies which have been 
documented by several authors 4 " 6 . The calculation of performance degradation due to 
icing is an equally important aspect of the icing analysis problem. Over the past several 


2 


years 2D calculations for iced airfoils have been performed using the interactive 
boundary layer method 7 and using Navier-Stokes codes 8 . Recently, Cebeci 9 has 
combined his interactive boundary layer method with LEWICE. This has resulted in a 
code which has the potential to accurately determine the flowfield and particle 
trajectories while also providing the performance degradation information during a 
single calculation. 

With these activities firmly underway, it is appropriate to consider the 
methodology for development of a three-dimensional ice accretion analysis capability. 
Some of the issues associated with a 3D icing calculation are the same as those for 2D. 
On the other hand, some issues only become evident when considering a full 3D 
calculation. These may include; variations in ice growth on both a large and small scale 
across a wing span, runback in the spanwise direction, non-uniformities in the impinging 
cloud along the span, and particle trajectories in complex geometric regions. 

The current effort is directed at taking existing computer codes and developing a 
strategy for performing a 3D ice accretion and performance degradation calculation for a 
swept wing configuration. This activity will serve as a demonstration calculation and as 
a focus for future development of the approach. The computer codes used in this effort 
consist of a 3D panel code, a 3D particle trajectory code, and the ice accretion 
subroutines of the 2D LEWICE code. The performance calculations will use a 3D 
Navier-Stokes code, ARC3D. Presently, the ice accretion and performance calculations 
are performed separately. In the future, these calculations will be performed in a single 
code thus requiring less user interaction than is currently necessary. 

This report addresses only the ice accretion aspects of the 3D problem. It 
documents the first attempt at comparison of 3D ice accretion calculations to measured 
ice shape tracings. The wing geometry chosen for this calculation is that of an MS-317 
airfoil with a 30 sweep angle, as shown in Figure 1. As indicated in the figure, the 
MS-317 profile was maintained in the direction parallel to the free stream velocity and 
not in the direction perpendicular to the leading edge. The MS-317 airfoil model 
represents a typical medium-speed wing section. The ice tracings were taken from a 
swept MS-317 airfoil model which was subjected to a range of icing condi tons in the 
NASA Lewis Icing Research Tunnel. The details of the test procedure are outlined in 
reference 10. 
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II. Code Description 


A 3D ice accretion calculation for a complete aircraft is a major step forward 
from current capabilities. The present activity is directed at extending current 2D 
capabilities to simple 3D geometries by assuming constant spanwise gradients. For this 
work, an infinite swept wing of constant geometry in the spanwise direction is used. 
Such a geometry produces similar streamlines along the span. As a result of this 
streamline similarity, it is assumed that there is no net transport of water across the 
control volume in the spanwise direction. The control volume analysis of the ice 
accretion can thus be considered two dimensional and the 2D LEWICE routines can be 
used. This also allows a single evaluation of the ice growth for any streamline to 
represent the ice accretion behavior for any location on the wing. Since the resulting ice 
shape geometry will not vary along the span, re-paneling of the iced wing for subsequent 
ice accretion calculations is simplified. 

An infinite swept wing also avoids consideration of ice growth near the tip of the 
wing or at the wing-body junction. At these locations, the flowfield is quite irregular 
and major changes to the approach currently being employed would be required. It is 
anticipated that future work in 3D ice accretion prediction would be directed at 
addressing these regions of the flowfield. Currently, the tools necessary for evaluation of 
these complex flowfields are being developed by Sankar 11 and Caruso 12 with funding 
from NASA Lewis. Additional work will be required to create a 3D control volume 
analysis for the ice accretion calculation. As these codes become available, they will be 
incorporated into the 3D ice accretion analysis methodology. 

II- 1. 3D Panel Code and Particle Trajectory Calculation 

The flowfield and particle trajectory calculations are extensions of the 2D 
analysis methods used in LEWICE. The three-dimensional flowfield calculation uses a 
Hess-Smith panel code 13 . The code can accomodate lifting and non-lifting geometries or 
combinations, thus allowing simulation of complex surfaces such as a complete aircraft. 
A Prandtl-Glauert correction is used to allow compressible flow calculations. The code 
can also handle leaking panels to simulate inlets or instrument orifices. The code 
provides velocity information at any point in the flowfield away from the surface. This 
type of code uses very little computer time compared to a viscous flow code and yet 
provides sufficient flowfield information for subsequent trajectory and energy balance 
calculations under most circumstances. 
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The trajectory calculation is based on the computer code developed by 
Norment 14 with one additional feature. The code solves a force balance equation on a 
single particle and then determines the new particle position using an Adams-type 
predictor-corrector algorithm developed by Krogh 15 . The added feature is the ability to 
calculate the local collection efficiency from the impacting particles. The code generates 
an array of impingement points for each region of interest. This is done by releasing an 
array of particles which impact the region of interest (Figure 2). The collection 
efficiency is the ratio of the particle flux at the target point to the free-stream particle 
flux. If a known group of particles is tracked from their release location to the target 
location, then the collection efficiency can be determined by the ratio of the target area 
to the release area. The collection efficiency can then be determined from the following 
relationship, 


/?(x c (i j)) = A o(bj)/A m (i,j) 


( 1 ) 


Once the collection efficiency is determined for the surface, an array of surface 
locations and corresponding /? values is created for streamlines on the surface. This 
information will be used in the ice accretion calculation to determine ice growth along 
the streamline paths. The streamlines are calculated using a 4th order Runge-Kutta 
integration scheme. The streamline is carried forward from the stagnation region for 
both the upper and lower surface at the region of interest. 

A linear interpolation scheme is used to determine the collection efficiency along 
the streamline from the matrix of /? values generated in the trajectory calculation. The 
array of /? values is searched to find the surface cell in which the streamline point 
resides (Figure 3). The /? value at the streamline point is then determined from a 
weighted average of the local collection efficiencies at the four surrounding cell points. 
The interpolation expression is, 


/?(x s ) ~ /?(x c (i+l,j+l))-u-v _|_ /?(x c (i,j+l))-u-(l-v) + /?(x c (i,j))-(l-u).(l-v) 

+ /?(x c (i+l,j))-(l-u)-v 


( 2 ) 


where, 



( 3 ) 
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and 


v = X. • Xj 


(4) 


represent weighting factors in the interpolation based on the geometry of the cell. 

II-2. Ice Accretion Calculation 

The ice accretion calculation is performed as a 2D strip analysis along 
streamlines calculated by the 3D panel code. The results of the previous calculation are 
written into a file containing the x,y,z coordinates of the streamline, the surface 
distance between points, the x,y,z components of the unit normal vector for the panel at 
the streamline coordinate locations, the velocities, and the local collection efficiencies, /?. 
The subroutines from the LEWICE code which determine the boundary layer values, 
the control volume energy balance and the resulting ice growth are used to determine 
the ice shape at a given spanwise location. A brief outline of the LEWICE calculation 
follows. Further details can be found in reference 3. 

The ice accretion calculation consists of a control volume mass and energy 
balance. The 2D strip analysis allows use of the assumption that there are no spanwise 
fluxes into the control volume, which is consistant with the infinite swept wing flowfield 
assumption. A depiction of the control volume for the process is shown in Figure 4. 
The As value is based on distance along the streamline. The surface area covered by 
the control volume consists of a parallelogram and is given by, 

A, = As c W (5) 

where As c is the component of the segment length in the chordwise direction and W is 
the width of the segment in the spanwise direction. A series of these control volumes 
are stacked along the streamline starting at the leading edge of the wing. 

The mass balance equation is, 

m c + m r . = m e + m r + m- (6) 

The impinging mass flux, m c , is determined from the /? distribution. The mass flux 
into the control volume, 1 %.^, is any water remaining in the liquid state after 
evaluation of an upstream control volume. This mass flux is called runback. The mass 
flux out of the control volume due to evaporation, m e , and to freezing, rh-, are 
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determined from the energy balance calculations. The mass flux out of the control 
volume, hi rout , is used as input for the downstream control volume. The freezing 
fraction is a term used to describe the fraction of incoming liquid that freezes within the 
control volume. It is given by the expression, 


/ = 


m- 


rii c + rii r ,n 


(7) 


The energy balance equation is, 

Qc d" qr in Qcond d" Qevap d" q conv d” Qr ou ^ d~ 


( 8 ) 


The values of q c and q r are obtained from mass flux values and the internal 
energy levels of the incoming water. Integral boundary layer expressions are used to 
determine the heat flux due to convection. The inviscid flowfield values are used as 
boundary values for the viscous calculation. The boundary layer calculation determines 
the location of transition from laminar to turbulent flow and then uses an appropriate 
value of the convection heat transfer coefficient. The heat transfer value is then 
obtained by use of Newton’s law of cooling, 


W = h A g (T, - T^ (9) 

The evaporative heat flux is also obtained from consideration of the inviscid flowfield 
conditions using the expression, 


1 evap 


= m # 


( C P... T. 


W, 8 


+ 


L„)as 


( 10 ) 


where the evaporative mass transfer rate, m e , is determined from local pressure, 
temperature, and humidity conditions. Further details on this calculation are found in 
reference 3. The heat flux due to conduction is a user specified value. This leaves two 
equations for the two unknowns, m r and m-. 

x out 1 

The mass of ice that freezes is used to determine the thickness of the resulting 
ice layer. Since the analysis is locally 2D, the ice thickness can be found from the 
relation, 


m-At 

d.„ =- 1 


Pice As c W 


( 11 ) 
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The thickness is then considered to be uniform over the entire panel for subsequent 
flowfield calculations. The new coordinates for the panel are obtained from the relation, 

x t - = x, + d, ce x,- (12) 

where x, is the coordinate of the center of the panel in the i-direction and x, is the 
i-component of the unit normal vector for the panel. 

As the ice thickness increases, there is the possibility that ice segments will 
intersect and thus this must be accounted for in the determination of the new geometry. 
Since this is a strip analysis, the ice thickness does not vary along the span at a given 
chordwise location. Therefore, the possibility of ice growth intersection is limited to the 
normal and chordwise directions. In that case, the line segments corresponding to the 
top of every other panel are examined for intersection. If intersection is determined to 
occur, then a new panel is formed with its center halfway between the two old panels. 
This requires determination of the coordinates of the new panel and renumbering of the 
panels. This information is then used in subsequent potential flow calculations. 

The ice shapes developed from the current calculations are based on a single time 
step. This was done to limit the number of particle trajectory calculations in favor of 
evaluating a larger number of icing conditions. The results, described in section IV, 
indicate that for short icing encounters the code compared well with experiment while 
for longer icing encounters time stepping may be required. In such cases, the output 
geometry file of the code is in the appropriate format for development of new surface 
panels. The updated geometry would then be used in a recalculation of the flowfield, 
particle trajectories, and ice growth. 

III. MS-317 Ice Shape Tracings 

The MS-317 airfoil used to generate ice shape tracings for a swept wing 
configuration is shown in figure 5. This model was created specifically to produce ice 
shapes for code validation purposes. As such, it was designed with a constant MS-317 
profile in the flowfield direction at every spanwise location, as opposed to a constant 
MS-317 profile perpendicular to the leading edge. This means that the chord of the 
airfoil is in the direction parallel to the tunnel walls and not perpendicular to the 
leading edge. The wing was also a full-span model in order to avoid the complications 
of tip effects on the ice growth patterns. The icing conditions tested, listed in Table 1, 
cover a range from glaze to rime conditions with velocities appropriate for a medium- 
speed airfoil. 
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Ice shape profiles were documented by cutting a slice in the accreted ice, parallel 
to the free stream flow direction, and tracing the ice profile onto a template. The cut in 
the ice shape was made at the mid-span location to ensure that end effects were 
eliminated and that cloud uniformity was maximized. The tracings were then digitized 
for use in plotting the profile and for comparison to predicted ice shapes. 

Since ice shapes were traced at only one spanwise location, the variability of the 
accreted ice in the spanwise direction is not documented quantitatively. Photographs of 
the ice accretions were taken however, as seen in the examples shown in figure 6, and 
they indicate a high degree of uniformity along the span in terms of ice mass and icing 
extent in the chordwise direction. There is however, considerable spanwise variation in 
the ice shape profiles. This type of variation in ice shapes has been documented in 
previous studies 16 and is demonstrated in figure 7. Also, the variation in chordwise 
extent of the ice is seen in the rime ice accretion photograph. This localized ice shape 
variability essentially precludes exact reproduction of an ice shape by any ice accretion 
calculation method. Instead, the code should be able to reproduce major features of the 
ice accretion, such as mass of the accreted ice as measured by the cross-sectional area, 
ice thickness at corresponding locations along the surface, presence or absence of horns, 
and average chordwise extent of the accretion. 

IV. Code Prediction/Tracing Comparisons 

The ice shapes were accreted at 2° and 8° angle of attack. The calculated 
pressure coefficient plots for those conditions are shown in figure 8. The plots show the 
large suction peak at the high angle of attack condition, which can alter trajectories of 
smaller droplets and tends to freeze out any water impinging on that surface due to the 
enhanced heat transfer. The relatively flat pressure distribution on the lower surface 
does not significantly alter droplet trajectories and allows greater runback on the lower 
surface. 

The calculated local collection efficiencies for the two angle of attack conditions 
are shown in figure 9. These results are for a droplet size of 20//m. The results indicate 
the reduced overall impingement at 8° and the shifting of the impingement pattern to 
the lower surface. The plots indicate that there is no impingement on the upper surface 
for the 8° angle of attack condition. Any calculated ice growth on the upper surface is 
therefore a result of freezing of runback water. These calculations are for a 2D cut 
along the flow direction. 
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The rime ice conditions correspond to the runs at 0°F (z.e. runs 8-11). The 
overall features of these ice shapes were reproduced very well by the calculations, as 
seen in figures 10-13. This is indicated by the impingement limit locations, which 
correspond to the ice limits for a rime ice growth, and by the cross-sectional area of the 
ice shape. Since water freezes on impact for rime ice cases, these results suggest that 
the flowfield, droplet trajectory, and ice growth routines are working well. The results 
also indicate closer agreement for the 6 minute results than for the 19.4 minute results. 
Since the code was run with only one time step for both cases, the changes in the 
flowfield due to the ice shape growth were not reflected in the collection efficiency or 
heat transfer calculations. Multiple time step calculations should be performed to try to 
improve the agreement between code and experiment. 

Mixed glaze and rime ice growths resulted from the conditions used in runs 1, 6, 
and 7. The comparisons are shown in figures 14-16. Once again, the agreement was 
much better for the shorter icing time. The a = 8° case, as shown in figure 16, indicates 
that the measured ice shape appeared to grow upstream as opposed to normal to the 
surface. The changes in collection efficiency distribution due to the development of the 
ice may be favoring growth in that direction. This interaction of changing particle 
impingement patterns with ice growth is not be accounted for by use of a single time 
step calculation. 

The runs at 28°F represent glaze ice accretions. These are characterized by clear 
ice, growth of horns, and runback water freezing aft of the impingement limits. Typical 
results are shown in figures 17-20. As seen, the calculations produce a shape that 
indicates a somewhat larger amount of runback than is evident from the actual ice 
growth. The calculation did not result in development of horns, which is due to the 
interaction of collection efficiency changes and runback. However, the overall growth 
patterns of the ice were reproduced quite well. 

The variability in ice shapes between two runs at the same condtions can be as 
great as the differences between measured and calculated shapes. This is due partly to 
uncertainty in the repeatability of set conditions in the actual cloud and partly to the 
spanwise variation of a given accretion. Figure 21 indicates this variability by 
comparing tracings at the same spanwise locations from two separate runs (i.e. runs 2 
and 13), nominally at the same conditions, to a calculated accretion. The measured 
shapes have a larger amount of total mass, most likely due to the increase in total 
collection efficiency for the actual accretion not accounted for by the single time step 
approach used in the calculation. On the basis of profile alone however, the calculated 
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shape more closely resembles the tracing from run 13 than that from run 2. This figure 
suggests that there may be little gained by trying to fine tune the ice accretion code 
beyond the point of determining the overall mass, whether it is rime or glaze, and the 
size and direction of glaze horns. The end product of the calculations, whether ice 
shape or performance degradation information, will also determine the degree of fidelity 
required from a given ice accretion calculation. 

The most influential parameter for glaze ice accretions, after the collection 
efficiency, is the roughness. The LEWICE roughness parameter establishes a roughness 
level used in the boundary layer calculations which is based on the liquid water content, 
velocity, and temperature of the cloud and the base surface roughness level of the wing. 
This roughness parameter was obtained from a correlation of these parameters with 
calculated ice shapes which compared well with measured ice shapes. This correlation 
was developed by Ruff and is described in the LEWICE user’s manual. 3 As the result 
of a series of calculations with different roughness levels, a roughness parameter a factor 
of four larger than the values suggested in the manual was found to produce better 
agreement with experiment and was used consistently throughout the course of this 
investigation. The need for a larger roughness value could be due to the spanwise flow 
and its effect on the heat transfer coefficient. Enhanced heat transfer due to this 
component of the flow would result in more rapid freezing and less runback than for 
two-dimensional airfoil models, from which the LEWICE correlations were derived. 
Thus, a new set of correlation equations may be required for ice accretion calculations 
on 3D geometries. This could prove a challenging task for more complicated geometries 
such as wing-body junctions. 

IV. Conclusions 

A 3D ice accretion code has been developed which calculates the 3D flowfield, 3D 
particle trajectories, and ice growth along 2D strips on the surface of interest. 
Comparisons to actual ice shape tracings indicate that the methods used in the 2D 
LEWICE code should be capable of being extended to evaluation of 3D ice growth. The 
single time step calculations used for this investigation were able to reproduce the 
overall characteristics of the ice shapes quite well. Extension to a time stepping 
proceedure, as in the 2D code, should be investigated to determine if the ice horn 
growth for glaze ice conditions can be more accurately reproduced. 

The approach employed in this investigation is applicable to infinite swept-wing 
geometries. This approach should also be adequate, considering all the uncertainties in 
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ice accretion calculations, in regions on the wing away from the root or tip, where 
spanwise gradients axe not severe. Rime ice conditions, where water freezes on impact 
and the determination of flow direction of runback water is not required, should also be 
capable of being predicted using this approach. However, in regions where the flowfield 
is very complex and under conditions where runback may occur, a different approach 
may be required to evaluate the development of glaze ice accretions. This could be 
quite significant for the evaluation of wing tips and outboard sections of rotor blades. 

The development of new heat transfer correlations may also be required if the 2D 
LEWICE methodology is to be extended to 3D glaze ice calculations. This would 
require the evaluation of ice accretion from several 3D wing models, in order to increase 
the confidence level in the code to that of 2D LEWICE. Additionally, some method of 
characterizing the spanwise variation of ice growth on a 3D geometry must be 
developed. This would allow more comparisons to code for a given icing experiment 
and indicate the degree of accuracy to which the code can be reasonably expected to 
conform. 
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Table I : Test conditions for MS-317 Ice Accretion Experiment 


Run 

T 

(T) 

V 

(mph) 

d 

(H 

LWC 

(g/m 3 ) 

t 

(sec) 

a 

( de s) 

(k/c ),- 

1 

15 

150 

20 

1.03 

1164 

2.0 

0.00192 

2 

28 

150 

20 

1.03 

1164 

2.0 

0.00931 

3 

28 

150 

20 

1.03 

390 

2.0 

0.00931 

4 

28 

150 

20 

1.03 

1164 

8.0 

0.00931 

5 

28 

150 

20 

1.03 

390 

8.0 

0.00931 

6 

15 

150 

20 

1.03 

1164 

8.0 

0.00192 

7 

15 

150 

20 

1.03 

390 

2.0 

0.00192 

8 

0 

150 

20 

1.03 

390 

2.0 

0.00127 

9 

0 

150 

20 

1.03 

390 

8.0 

0.00127 

10 

0 

150 

20 

1.03 

1164 

8.0 

0.00127 

11 

0 

150 

20 

1.03 

1164 

2.0 

0.00127 

12 

28 

150 
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SECTION a - a 

FIGURE 1. - MS-317 SWEPT WING PROFILE. 


+ 


Impingement point array C x m < I , J ) ] 


starting point array t x Q < I . J > ] 



FIGURE 2. - ILLUSIRAl ION OF THE STARTING POINT ARRAY lx 0 (I.J)J AND IMPINGEMENT POINT ARRAY (X m (I.J)l 
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Impingement point grid, x^Ci.j) 
--^-impingement centroid grid, x_ < i , j ) 



x c (i. 



x c (i+l,j+l) 


|3(Q) = P(x (i + i ,j + l))*u*v + 3<x (1 ,j+l))*u*(l-v) + P(x (i ,j))*(l-u)*(l-v) 

+ (3( x c ( i + 1 ,j)) • ( 1 -u ) *v 

FIGURE 3. - ILLUSTRATION OF THE CENTROID IMPINGEMENT GRID C S c ( I , J) 1 AND 0 INTERPOl ATION. 
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FIGURE 5. - MS-317 SWEPT WING MODEL USED TO OBTAIN ICE SHAPE 
TRACINGS. 



(a) RIME ICE AT 0 °F AND a = 2°, RUN NO. 8B. <b) GLAZE ICE AT 28 °F AND a = 2°, RUN NO. 13. 

FIGURE 6. - TYPICAL ICE ACCRETION ON MS-317 MODEL. TRACING OF ICE SHAPE WAS MADE FROM CUT PARALLEL TO THE FLOW DIRECTION. 
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Section A 



Section C z 


Section B 



Section D 


Section C, 



Section E 


(a) Ice shape tracings for sections a-e. 



icing band 



(b) Composite of ice shape tracings showing icing band. 

FIGURE 7. - SPANWISE VARIATION IN ICE SHAPES FROM A 2D ICE ACCRETION. TRACINGS TAKEN AT SEVERAL LOCATIONS ALONG THE 
SPAN OF A BOEING 737 MODEL TESTED IN THE NASA LEWIS IRT. (TAKEN FROM REF. 15.) 


19 
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FIGURE 8. - CALCULATED PRESSURE COEFFICIENT DISTRIBUTIONS 
FOR MS-317 WING. 



FIGURE 9. - CALCULATED COLLECTION EFFICIENCY DISTRIBUTIONS 
FOR MS-317 WING. ANGLE OF ATTACK, a, IS 2° AND 8°. 


2 — 



-2 — 




4 1 1 1 1 1 

-2 0 2 468 

X, IN. 

FIGURE 10. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 8.) 



X, IN. 


FIGURE 11. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION 
(RUN NO. 11). 



FIGURE 12. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION 
(RUN NO. 9). 
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EXPERIMENTAL 

CALCULATION 



EXPERIMENTAL 

CALCULATION 



■4 1 1 1 1 1 

-2 0 2 4 6 8 

X, IN. 

FIGURE 13. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 10). 
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FIGURE 16. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 6). 
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-2 0 2 4 6 8 

X, IN. 

FIGURE 14. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 7). 
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FIGURE 17. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 3) 
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FIGURE 15. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 1 ) . 
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FIGURE 18. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 2). 
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FIGURE 19. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO. 5). 



X , IN. 


FIGURE 20. - COMPARISON OF EXPERIMENTAL ICE SHAPE 
TRACING WITH 3D CALCULATION FOR MS-317 WING SECTION. 
(RUN NO.M). 



FIGURE 21. - COMPARISON OF EXPERIMENTAL ICE SHAPE TRAC- 
ING WITH 3D CALCULATION FOR MS-317 WING SECTION. VAR- 
IABILITY OF THE MEASURED ICE SHAPE IS INDICATED BY COM- 
PARISON TO TRACINGS FROM TWO RUNS AT THE SAME CONDITIONS 
(RUNS NO. 2 AND NO. 13). 
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